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Abstract- The performance of coherent Doppler lidar is determined by the statistics of the coherent 
Doppler signal. The derivation and calculation of the covariance of the Doppler lidar signal is presented 
for random atmospheric wind fields with wind shear. The random component is described by a Kolmo- 
gorov turbulence spectrum. The signal parameters are clarified for a general coherent Doppler lidar sys- 
tem. There are two distinct physical regimes: one where the transmitted pulse determines the signal 
statistics and the other where the wind field dominates the signal statistics. The Doppler shift of the signal 
is identified in terms of the wind field and system parameters. 


I. Introduction 

Coherent Doppler lidar has become a promising new instrument for observations of the atmosphere. 
The first systems were based on the CO 2 laser 1-4 . More recently, solid state lasers have been successful 
integrated in Doppler lidar designs 5-7 . The targets for Dopper lidar are small atmospheric aerosol parti- 
cles. Doppler lidar data can be the return from a single pulse or a collection of many pulses along the 
same line of sight. The advantages of multiple pulse has been discussed by Rye and Hardesty 8 , and 
Frehlich 9 . Coherent Doppler lidar data is used to estimate the radial velocity of the targets using various 
mean-frequency estimators 8-14 . 

The actual measurement quantity is related to the statistics of the random wind field over the sens- 
ing volume of the pulse. The effects of the pulse volume averaging for Doppler radar has been investi- 
gated by Doviak and Zmic 15 . The estimated Doppler frequency is defined in terms of the spectrum of the 
signal. Doppler radar data is generated by transmitting many pulses and sampling the returns at fixed 
sampling interval T s [sec], which corresponds to sampling the same spatial volume defined by the pulse 
at equally spaced time intervals. For this case, the statistics of the data are determined by the time evolu- 
tion of the sensing volume of the pulse, and the spectrum is a well defined quantity. For coherent 
Doppler lidar, the data is generated from a single pulse as it travels through the atmosphere. The sensing 
volume of the measured Doppler frequency is determined by the distance the pulse travels in the finite 
observation time of each estimate T=MTs, where M is the total number of data points per estimate. The 
spectrum is not well defined and the Doppler frequency is related to the statistics of the wind field over 
the sensing volume of the single pulse. The estimated parameters are then related to the wind fields and 
lidar system parameters in a different fashion than for the Doppler radar. This connection will be derived 
to clarify the meaning of coherent Doppler lidar measurements and pave the way for more reliable esti- 
mators of the signal parameters. 

The statistical description of coherent Doppler lidar data from aerosol targets has been shown to be 
well represented as a zero mean Gaussian random process. It is a common practice to follow the tradition 
of Doppler radar and convert the real signal to the complex analytical signal using a quadrature 
receiver 5-17 . This generates complex data with desirable statistical properties that simplify the calcula- 
tion of the Cramer-Rao Bound 14,15,18 (the ideal performance of estimators of the signal parameters) and 
the development of efficient estimators of the signal parameters. The statistical description of the data is 
required for these analysis. Since the data is a Gaussian random process, the signal covariance provides a 
complete statistical description of the data. The signal covariance is derived for a general coherent 
Doppler lidar including random wind fields and wind shear. The derivation follows the analyses of 
Frehlich and Kavaya 19 for the performance of coherent laser radar for general atmospheric refractive tur- 
bulence. 


2. Coherent Doppler Lidar Signal 

A coherent Doppler lidar can be represented by the transmitter aperture function WY(if), the 
transmitter laser field incident on the transmitter aperture E L (Ji,0,t), the receiver aperture function W R (V), 
and the Local Oscillator (LO) field incident on the receiver aperture Eiq (V,0). The transmitter laser field 
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with frequency ft [Hz ] is given by 

E L QS,0,t)=A L (t)e L (ji,0)exp(2Kifit) ( 1 ) 

where /[s ] denotes time, e L (tf ,0) is the nonnalized spatial field which is assumed independent of time t 
(the spatial mode of the laser pulse is independent of time) and 

/ 

0 

is the temporal behavior of the complex amplitude of the pulse where f*z,(0[W] is the pulse power and 
/c(OI7fc] is the frequency chirp of the transmitted pulse. The LO field in the receiver plane with fre- 
quency fto [Hz ] is given by 

Elo m = P$e w (V,0)exp(27u/ LO r) (3) 

where Pto [W ] is the LO power and 

oo 

J I e L to WO) | dlf=l (4) 

The coherent Doppler lidar typically operates under the ideal condition of quantum limited detection, i.e., 
the signal noise is dominated by the shot noise of the LO field. We will assume that the detector quantum 
efficiency x\q is uniform over the surface of the detector and that all the LO power is collected by the 
detector. 

For Doppler lidar measurements of the atmosphere, the backscattered field is a zero-mean Gaussian 
random process because it is the superposition of all the backscattered fields from many random aerosols 
is the superposition of the backscattered field from many aerosol particles. The doppler lidar signal is also 
a zero-mean Gaussian random process. It is convenient to normalize the signal by the standard deviation 
of the noise signal and employ the complex representation (Appendix A) 

Z(t)=S(t )+N (t )=Re[z (f )] (5) 

where S (r) is the normalized Doppler lidar IF current signal from the backscattered field and N(t) is the 
nonnalized additive detector noise. The statistical description of a Gaussian random process is given by 
the covariance function 


B z (t i ,t 2 )= <2 (t i)Z (t 2 )> = Re[/? z ( t j ,t 2 )] (6) 

**('i.'2)=y<z(fi)z*(f2» (7) 

<z{jt x )z{t^>=0 (8) 

where the ensemble average <> is performed over the appropriate random process. With this normaliza- 
tion 


fiz(fl>*2) = S(tl - t2) + 5.s(fl,f2)> Rz(fl.f2)=5( f l - *2H/?z(fl.f2) (9) 

where B s and /?.,(/ ^ 2 ) are the covariance of the normalized Doppler lidar signals, 

B s (t,t) = R s (t,ty=SNR(t) (10) 


and SNR is the coherent lidar signal-to-noise ratio with respect to the ensemble average. The Doppler 
lidar signal is typically uncorrelated from shot to shot because the aerosol targets move more than a 
wavelength between shots, which produces random phases for the backscattered signal. The ensemble 
average over the aerosol particles is the most common ensemble average and produces (Appendix A) 


J Al (' i-2|MZ(^-2^-)^ 2 (?.z W-z )h2(fT,z)| 2 exp|2ra x[AF--^ v, (p,z )] jdpdz (11) 


where v r (j t,z) is the radial component of the velocity of the aerosol particles at coordinate (ff.z), 
A F=ft~fbo is the frequency difference between the transmitted field and the LO field, and 
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Q(pj) = \\ jw r (i t)Wn (?)e L (!f,0 )e* w <?,0)G (phlf^ X? )</I?<f? (12) 

where G (pjt,z ) is the Green’s function for propagating the field from the transmitter plane to (J f,z). For 
Doppler lidar, the complex signal z (f ) is also an analytic signal because Eq. (8) is satisfied (see Appendix 
A). This provides simplifications for data analysis and estimation algorithms. The function R (ft, * 2 ) is a 
hermitian function, i.e., 

/?(f,,f 2 ) = /?*(f 2 ,f 1 ) (13) 


With a high pulse rate Doppler lidar, many shots can be transmitted during the temporal scale of the 
random Green’s function G (Jf;I f,z), which is usually determined by the atmospheric refractive turbulence. 
This permits investigations into the statistics of refractive turbulence. The ensemble average over refrac- 
tive turbulence produce 


R s(hJ 2 >=-r^i j a l fr !-2 z !c )A* l ( t 2 -lzic )0 (j?,z )exp| 


h\B 


2ni t[AF -2v r (Jf.z )/X] \dpdz 


(14) 


0(p,z) = K 2 (p,z)^)c(p,z), 


(15) 


where 


c(ff,z)= |G(p.z)l 2 . ( 16 ) 

is the coherent responsivity density 19 of the coherent Doppler lidar. For many applications, the effects of 
refractive turbulence are negligible and 

G(iJ;tf,z) = G/(jpW,z) = ^exp[--g-(jMi) 2 ] (17) 


where G f (jj;lf,z ) is the free space Green’s function and k=2n/‘k[m~*] is the wavenumber of the laser 
field. The fluctuations of the wind field ^(p,z), atmospheric extinction K (p,z ), and backscatter 
coefficient fKftz ), are the remaining random quantities. 

A typical Doppler lidar has a narrow beam compared to the dimensions of the range gate. Then, the 
dependence on the transverse coordinate 1? is small and 


R s (fi,f 2 >=-^-ltf(z )A l (f !-2z !c )A' L (t 2 -2zIc )expj 2nix[AF-2v r (0,z)/X] 


hvB 


(18) 


where 


is the system gain and 


H (z )=K 2 (0j )|3(0,z )C (z ) 


(19) 


C(z)=lc(p,z)dp ( 20 ) 

is the coherent responsivity 19 of the coherent Doppler lidar. 

The center of a range gate for Doppler lidar measurements is defined by 

zo=cf(/2 (21) 

and the total observation time corresponding to this range gate is (tg-T /2,t 0 +T 12). If Pflfcz ), K (| 1,z ), 
v r (|f,z), and c(p,z) are constant over the sensing volume of this range gate and if there is no frequency 
chirp /c(f). we have for all (f in the observation interval 

R s(tx,td=R s W = (z o)exp(27u' xF ) J [P L (t)P L (t-X)] m dt (22) 


where 
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F=bF-2v r (Oj 0 )l\ (23) 

is the mean frequency of the Doppler lidar signal, v=t\-t 2 and the limits of integration have been 
extended to infinity. The signals x{t) and Z(t) are stationary over the observation interval 
(r 0-T /2,t 0 +T 12). For a Gaussian pulse with 

P L (t>P<fixp(-t 2 /o 2 ) (24) 

R s (t) = SNR (z 0 )exp(2jc/F ~^/x P 2 ) (25) 

where 


X\q cUt 

SNR(z^-£——H(z 0 ) 


2hvB 


is the SNR for a pulse centered at range z 0 . 


U L =jP L (t)dt 


(26) 


(27) 


is the pulse energy and x P =2a is the decorrelation time of the lidar signal for uniform wind fields. 
The signal spectrum is defined as 


(/ ) = \ R i (x)exp(-2jtt' xf )dx 


and 


wTiere 


*,(/) = 


^-K\zo)^z 0 )C(z 0 )M(f-F) 


(28) 


(29) 


M(f)= | jP L (t) v hxp(-2mft)dtT 


(30) 


is the spectrum of the transmitted pulse. The signal spectrum is also the same as the signal spectrum from 
a rigid or diffuse hard target 

3. Effects of Wi nd Sh ear 

Doppler lidar data permits estimation of the radial component of the wind field v r (p,z ) with a single 
shot For short pulse duration, the random fluctuations of v r flf.z ) can be approximated by the first two 
terms of it’s Taylor series expansion in z , i.e., 

(31) 


V r (p,Z ) = v r (Jf,z 0 )+g (J*,z 0 Xz— Zo) 


where g df,z 0 )[5 ! ] is the wind shear of the radial component of the velocity at flJ,z 0 ), i.e., 


g(p.z) = 


jVrOM 

dz 


(32) 


The range dependent system gain 0(p,z) is approximated by the first two terms of it’s Taylor series 
expansion in z, i.e., 

O (i?,z ) = O Of,zo)[l+& (if,z 0 )(z-zo)] (33) 

where - — 

_ i aoqfc) 

0(tf,z) dz 

When these approximations are valid and for a Gaussian pulse [see Eq. (24)] 


b(p,z) = 


(34) 
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Rsih,t ^ = ~^t J O^z 0 )[\+2b(^a Q W{\U<y-t^<5-i^—)] 


exrt 


2ni x[A F -2v r (J f,z )/X+2(.t (r M'tws o)]-x*/x F 2 ^i 2 rt ws 2 W]t 


where \i^(ti+t?)/2. 


tws = 2 


g(p,z 0 )nAr 


! 2xws 
■% 



(35) 


(36) 


(37) 


is the decorrelation time due to wind shear and 

Ar=col2 

is the length of the pulse in range or the range resolution. For a Doppler lidar with a narrow beam com- 
pared to the dimensions of the range gate 


R s (t u t£=SNR (z 0 )[ l+2'y(n/a-r(/o-/^— )] 
expjz t[2tcF -2(p-r 0 )/(xws o)]-^n P 2 ^/x ws 2 


(38) 


where 

y=b (0,z 0 )Ar =[0 (0,zo+Ar )-0 (0,z 0 )]/O (0,z o) (39) 

is the fractional change of the gain O (0,z ) over the range resolution Ar at zo, if the linear approximation 
of Eqs. (31) and (33) are valid. 

The rate of decay of the covariance is dominated by the wind shear when 

2rcc a 2 # (0,z o )/A» 1 (40) 

This occurs when the extra Doppler shift due to the wind shear at a separation of the range resolution Ar 
is larger than the spectral width of the signal due to the pulse profile only [see Eq. (29)] 

The effects of wind shear with constant SNR [H(z)] over the observation interval are shown in Fig. 
1. The wind shear produces a rapid decorrelation of the lidar signal, which appears as a damped oscillat- 
ing signal as the observation points ti and t 2 deviate from the center of the observation interval t 0 . The 
effects of a gradient in the SNR [//(z)] is shown in Fig. 2. The gradient produces an amplitude change in 
addition to the effects of wind shear. 


3- Effects of Wind Turbulence 


Atmospheric wind fields are characterized by random fluctuations in space and time. For many 
cases, the statistics of the wind fields are well characterized by universal functions of basic parameters. 
For Doppler radar the effects of the random fluctuations of the wind field over the sensing volume has 
been discussed by Doviak and Zmic 15 . For Doppler lidar, the situation is complicated because estimates 
of the velocity are computed for a single shot. Ensemble averages over many shots was considered by 
Chumside and Yura 21 . The effects of spatial variations of backscatter coefficient over the sensing volume 
of the pulse were investigated by Rye using simulations 22 . 

The covariance of the Doppler lidar signal for a single shot is given by Eq. (14), which can be writ- 


ten as 


R s (t \J2 >=y^r J J ^(r t -2z/c W(f 2 - 2 z/c) 0 (ir,z)exp| 


h\B 


2iUT[F-2Av r (ff,z)/X]} 


(41) 


where 

Av r (jf,z ,zq) = v r (jfc )-v r (0,z 0 ) ( 42 ) 
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is the difference of the radial component of the velocity compared to v r (0,z 0 ). The mean frequency / is a 
random variable because the wind field is random. The estimation error for estimates of the mean fre- 
quency F using data from a single shot depends on the statistical properties of Av r (j/z ,z 0 ) and the func- 
tion H (jf.z ). For well behaved atmospheric conditions, the velocity difference Av r (fJ,z,z 0 ) is approxi- 
mately a zero mean Gaussian random variable. Performing the ensemble average over this random vari- 
able produces 

oo (JO 

i.'2>=-^r J j a l i~2z fc )A* L {t 2 -2z Ic )H (j/z )exp[27ri tF-2t 2 D, (jj,z ,z 0 )/X 2 ]dz (43) 

where 


D r (I t,z ,z 0 ) = <Av r (| t,z ,zo) 2 > 


(44) 


is the structure function of the fluctuations of the radii component of the velocity field. When the statis- 
tics of the velocity fields over the sensing volume of the pulse are isotropic 20 

D r (jt,z ^o) = C v e 2/3 (z 0 )r 2/3 [(z-z 0 ) 2 44p 2 /3]/r 2 (45) 

where C v is the Kolmogorov constant, e[m 2 j -2 ] is the energy dissipation rate, and r^z-z o) 2 +p 2 . Fora 
Guassian transmitted pulse 

TV>f//cAr 2 „ 

Rs (fi,f 2 >= — 1/2 exp(2 iu xF-^IXp 2 ) 

oo oo r ^ 

J J H (Arif,Ar£+zo)exp -[(|i-r 0 )/a-^] 2 -8T 2 AK 2/3 (^ 2 44^ 2 /3)/(KT/>) 2 Lf^ii^ (46) 

where k 2 =£ 2 +< 7 2 and 

A=o 2 C v (eAr ) 2J3 /X 2 =o 2 D r (O.zo+Ar j 0 yk 2 (47) 

represents the contribution from wind turbulence. 

For a coherent Doppler lidar with a Gaussian transmitted spatial profile 19 

H (Ar <f,z 0 ) = exp[-Ar 2 q 2 /o£ (z a)] (48) 

where C (zo) is the coherent responsivity at range z o and Or ( z o) is the 1/e intensity radius of the transmit- 
ted beam at range z 0 . For many cases of interest, C (z 0 ) and o B are constant with range over the sensing 
volume of the pulse and 

2SNR(td)Ar 2 . , 

R s (tutz)= 7rT~i exp(27« XF -Xr/Xp 2 ) 


J J ex H 


rc ,/2 c£ 

-[(^-r 0 )/o-Q 2 -8x 2 AK 2/3 (C 2 +4< ? 2 / 3 )/( kx p ) 2 fd«jtf C 


For a narrow beam, the integration overscan be performed and 
SNR (z 0 ) 


R s 0 1»#2)= — — exp( 27 rixF^r 2 /T^ 2 ) J expj 


[(ft-f o)/<*-C] 2 -8T 2 Ai; 2/3 /T / , 2 \d C, 


(49) 


(50) 


The effects of wind turbulence and the finite transverse dimensions of the beam is shown in Fig. 3. 
The decrease in the correlation scale due to turbulence is more pronounced for a wide beam (a B >A r) 
because of the increase in the Doppler shifts from the random fluctuations of the wind in the transverse 
dimensions. 

The effects of wind turbulence for the common case of a narrow transmitted beam compared to the 
range resolution (c B <A r) is shown if Fig. 4. 
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5. Summary 

The signal covariance for coherent Doppler lidar data is derived including a linear wind shear and 
linear system gain over the sensing volume of the pulse. Results are also presented for the case of isotro- 
pic and homogeneous wind turbulence. For both these case, the signal statistics are not stationary over the 
observation interval. The use of estimators for the system parameters that assume stationary (spectral 
domain estimators) may not perform efficiently in these regimes. Improved performance is possible 
using the improved description of the data contained in the signal covariance. 
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APPENDIX A: Coherent Doppler Lidar Signals 

The total coherent or heterodyne Doppler lidar signal iV(f)[A] is the sum of the Doppler IF signal 
current /5(f) [A] and the detector noise current (v(0[A ] which is conveniently written in complex nota- 
tion as i T (t)=Re [T(f )], i N ( t)=Re [£!(#)], where 

Y(t)=U(ty+^^-^r\Q&)E s &£ ,t )Eio )exp(i Acm +i e s )d& (Al) 

Go is the dimensionless amplifier gain, e = 1.602 x \QT^[C I electron] is the electronic charge, 
h =6.626 x lCr 34 [/s] is Planck’s constant, \[Hz] is the optical frequency, r\ Q (^[electrons /photon ] is 
the detector quantum efficiency function on the detector surface, E s ,t ) is the backscattered field on 
the detector surface, Euo ) is the LO field on the detector surface, 85 [rad] is the random phase of the 
backscattered field compared to the LO field, &<n[radlsec], is the angular frequency difference between 
the transmitted and LO field, [m ] is the transverse coordinate on the detector surface and d\t denotes 
two-dimensional integration over the detector surface. For ideal optical systems and for a infinite uniform 
detector, the fields on the detector can be represented in terms of the fields in the plane of the receiver 
optics, i.e., 


Y(t)=U (t) + J E s <?,0,t )Eto (?,0) W* (V)exp(j A(M +i 0 S )d? (A2) 

where E s (?,0,f) and E L0 <V,0) are the backscattered field and LO field in the receiver plane. 

If the signal noise is dominated by the shot noise generated by the local oscillator field (quantum 
limited detection), the average noise power is 

<i$(t)>=j<U(t)U\t)> =2G D 2 e 2 BT\ Q P LO l(hv) (A3) 

where B [Hz ] is the noise bandwidth. It is convenient to normalize the signals by the rms of the noise 
current, i.e., 

Z (t y=Y(t)/<rf(t )> 1/2 = S(tyt-N (t ). (A4) 

The backscattered field from a point scatterer at coordinate (ftz) with scattering cross section c s is 

00 

E s (V,0,r ) = Xa} a K (Jf.z ) j E l (tf,0,r-2z /c )W T (fl)G (ff.Tf.z )G (jfcVj )exp[i 8(ff,z )-4iutv r (p,z )/X]dI? (A5) 

where 0Qp,z ) is the random phase of the backscattered field. The total backscattered field is the sum of 
all the backscattered fields from the aerosol particles illuminated by the transmitted field. Substituting the 
total backscattered field into Eqs. (A2), (A4), (7), (1), performing the ensemble average over the random 
phases 8flf,z) using the identity <exp[/ (©(pf.z 1 >+8(p 2 ,z 2)]>=0. and converting the summation over 
aerosol particles to integration 19 produces Eq. (1 1). The same procedure produces Eq. (8). 
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Effects of wind shear on normalized signal covariance [Eqs. (6), (9), (38), f =0] as a function of 

x/Zp with constant SNR over the observation interval and ( __) |i=ro, no wind shear, (....) 

p.— t(f=G and Zp=^iws'’ ( — ) ft — f(p2<j and Xp = 1 ^ 5 » (•■■■) ft— ro=lo and Xp=2xws\ and .. -) 
u — t o=2<j and Xp =2 Xws ■ 
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Figure 2 


Fig.2 Effects of wind shear and gradients in SNR over the observation interval on the normalized sig- 
nal covariance [Eqs. (6), (9), (38), f=0] as a function of xl x P compared with constant SNR over 

the observation interval and no wind shear ( ). The curves represent (....) x P =Xws , 

and y=0.2; ( ) p-fo=o, x P =x W s - and Y=0- 4 : (-•-•) \i-ta=2o, x P =2x ws , and y=0; (.. - .. -) 

|i-f(j=2o, x P =2xws • and Y=0-2; and |x-r o=2o, x P =2x ws , and y=0.4. 




Figure 3 


Fig.3 Effects of the finite transverse dimensions of a lidar beam and wind turbulence on the normalized 
signal covariance [Eqs. (6), (9), (49), ^=0] as a function of x/Tp with constant SNR over the 

observation interval and case a) jx=f o with (_ ) no wind turbulence; (...) A=2, and c B «A r ; (- -) 

A=2, and o B =Ar/2; (.-.-) A=2, and o B =Ar; (--.--.) A=2, and g b = 2A r ; and case b) the same 
as case a) but with ji-/(y=a. 
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Figure 4 


Fig.4 Effects of wind turbulence on the normalized signal covariance [Eqs. (6), (9), (50)*F=0 ] as a 
function of x/x P with constant SNR over the observation interval and narrow transmitted beam 

compared to the range resolution (o fi «Ar) for case a) |X=r 0 with ( ) no wind turbulence; (...) 

A=2; (- -) A=4; and case b) is the same as case a) but with |Wo=o. 
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